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ABSTRACT 


For  the  purpose  of  determining  a  general  correlation  for  forced 
convective  heat  transfer  coefficients  for  circular  cylinders  in  cross-flow 
under  conditions  pertaining  to  cooled-filir.  application  for  measurements 
in  flakes  and  hypersonic  wakes,  heat  transfer  measurements  were  obtained 
under  simulated  conditions  created  in  a  plasma-jet.  Measurements  were 
made  by  means  of  constant  temperature,  internally  cooled,  cylindrical 
film  sensors  in  the  Reynolds  number  range  5-75  for  various  temperature 
loadings,  the  maximum  being  T(„,  /'fg  =  4,  and  for  flows  composed  of  He, 

N2  and  CO2.  The  plasma-jet  conditions  were  maintained  such  that  at  the 
points  of  measurements,  i.e.  at  the  potential  core  of  the  jet,  ionization 
was  negligible,  recombination  was  complete  and  the  turbulence  intensity 
was  low.  3y  a  statistical  analysis,  the  data  obtained  for  Rem  »  5  to  40 
were  correlated  by 


.  .  •  T5  /c 

NV  *  -2068  +  *4966  (Re*) 

All  the  fluid  property  values  in  the  dimensionless  parameters  were  evaluated 
at  the  arithmetic  mean  temperature.  A  kinematic  viscosity  ratio  instead 
of  the  more  usual  temperature  ratio  was  considered  in  the  temperature 
loading  factor  to  account  for  the  different  dependency  of  the  transport 
properties  on  temperature  for  different  gas  species. 

Though  no  previous  correlations  were  obtained  for  conditions  con¬ 
sidered  here,  the  direction  of  the  temperature  loading  effect  and  the 
departure  from  King T 5  law  appear  to  be  consistent  with  scxne  of  the  existing 
correlations. 
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Greek  Symbols 
00 

>> 

§ 


NOMENCLATURE 


Constant 

Surface  area  of  the  cylinder  (sensor) 

Constant, 

Constant 

Specific  heat 

Diameter,  sensor 

Mass  velocity 

Heat  transfer  coefficient 

Thermal  conductivity 

Constants 

Nusselt  number  hd 
K 

Frandtl  number  c_  u 

ST 

Heat  transferred 
Resistance,  eiec+rical 
Reynolds  number  Vdp  ,  Vd 

r-  >> 

Temperature 

Fluctuating  component  of  velocity  parallel  to 
mean  flow  velocity 

Mean  flow  velocity  normal  to  cylinder  (sensor) 

Temperature  coefficient  of  resistance 
Absolute  viscosity 
Kinematic  viscosity 


Mass  Density 
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NOMENCLATURE  (Confd) 


Subscripts 


m  Arithmetic  mean  and  when  referring  to  dimensionless 

groups  signifies  that  the  fluid  property  values  are 
evaluated  at  the  arithmetic  mean  temperature,  i.e. 
at  Ts  +  T« 

2 

s  Surface,  sensor 


oO 


Free  stream 
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1.0  INTRODUCTION 

1.1  Direct  measurements  of  dynamic  characteristics  in  high  temperature 
fields  as  exist  in  laminar  and  turbulent  flames,  hypersonic  wakes,  rocket 
nozzle  flows  demand  a  rapid  response  instrument  able  to  withstand  high 
temperature  environments.  An  extension  of  constant-temperature  hot-wire 
anemometry  to  constant-temperature,  internally  cooled,  thin-fiim  anemom- 
etry  towards  this  end  is  considered.  (Compared  to  hot-wire  sensors, 
cooling  of  the  cylindrical  film  sensor  not  only  presents  melting  of  the 
sensing  element  but  also  allows  interpretable  measurements  to  be  obtained 
when  the  environment  is  at  a  higher  temperature  than  the  sensing  film,, 
without  sacrificing  response) .  The  principle  of  operation  of  the  cooled- 
film  sensor  as  compared  to  that  of  the  hot-wire  or  hot-film  sensors  is 
shown  in  Fig.  1. 

For  cooled-film  operation  when  the  environment  temperature 
is  greater  than  the  sensor  temperature,  a  heat  balance  applied  tc  the 
sensing  film  gives  (See  Fig.  1): 


=  Coolant  ~  Supplied  .  (1A) 

In  the  above  relation,  Q~.jpp-  j™  is  the  heat  equivalent  of 
the  electrical  power  supplied  by  an  anki&gnetwork  to  maintain  the 
constant  film  temperature  Ts.  Qgjjy  is  the  heat  absorbed  by  the  film 
from  the  environment  and  OcoOLANT  is  the  heat  removed  by  the  coolant. 

From  equation  1,  the  parameter  of  interest  Qenv  i£  obtained  by  measuring 
^COOLANT  SUPPLIED'  As  long  as  the  coolant  flow  rate  and  inlet 
temperature  remain  constant*  heat  absorbed  by  the  coolant  (QcOOLANT)  from 
the  constant-temperature  film  will  remain  constant  and  is  independent  of 
the  environmental  conditions.  The  actual  value  of  this  OcoOLANT  i® 
obtained  from  a  measure  of  the  QsUPPLIED  an  environment  where  Qgjjy  0. 
OsUPPMED’  a^  any  time,  is  obtained  by  measuring  the  current  in  the 
network  required  to  maintain  Ts  constant.  In  effect  then,  under  operating 
conditions. 


qenv  =  Supplied  (when  o^y-^o)  '  Supplied  (ib) 

The  Qgfjy  measured  in  this  way  then  may  be  interpreted  in  terms  of  the 
parameters  of  interest  of  the  flow  field  by  means  of  an  accurate  heat 
transfer  relation.  The  purpose  of  the  present  experimental  investigation 
was  to  cbtain  this  heat  transfer  relation  under  conditions  of  such  variables 
ranges  that  are  likely  to  be  encountered  in  tne  application  of  cooled-film 
sensors  to  laminar  and  turbulent  flames  and  hypersonic  wakes. 

1.2  Since  Qgjy  is  predominantly  due  to  forced  convection,  then 
the  forced-convective  heat  transfer  relation  for  single  cylinders  in  in¬ 
compressible  cross-flow  is  required.  Further,  th*  relation  must  be  valid 
for  tne  following  cases: 

a)  When  tne  transfer  of  heat  is  from  the  environment  t.o  the 
cylinder. 
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b)  When  the  temperature  difference  between  the  cylinder  and 
the  environment  is  large,  and 

c)  When  the  environment  i3  composed  of  flows  of  different 
gases  and  gas  mixtures. 

1.3  For  forced-convective  heat  transfer  between  single  cylinders 
and  incompressible  cross-flows  under  near-isothermal  conditions  (or  on 
the  basis  of  the  assumption  that  the  fluid  properties  are  independent  of 
temperature),  a  dimensional  analysis  of  the  pertinent  parameters  gives 
the  following  relation; 


(hd)  -  f(M  ,  )  (2) 

k  V  k; 

However,  for  large  temperature  differences  existing  between 
the  heat  transfer  surface  and  the  bulk  of  the  fluid,  when  a  considerable 
gradient  of  fluid  properties  exists  in  the  thermal  boundary  layer,  the 
above  dimensionless  groups  cannot  be  uniquely  defined.  For  the  purpose 
of  correlation,  this  so-called  temperature  loading  effect  must  then  be 
taken  into  account  by  evaluating  the  fluid  properties  at  suitable  mean 
temperatures  -  arithmetic  mean,  weighted  arithmetic  mean..  logarithmic 
mean,  geometric  mean,  or  even  the  individual  fluid  properties  at  different 
temperatures.  Since  the  temperature  loading  effect  is  empirical,  the 
choice  depends  solely  on  the  success  arid  convenience  in  the  evaluation 
and  application  of  the  correlation. 

However,  it  is  to  be  realized  that  for  non-isothermal  cases, 
a  consideration  of  the  dimensionless  groups  evolved  for  near-isothermal 
cases  alone  does  not  suffice  for  the  condition  of  dynamic  similarity. 

The  condition  of  dynamic  similarity  is  that  the  same  ratios  must  exist 
between  the  significant  fluid  properties  at  geometrically  equivalent 
points.  Since  the  gradient  of  properties,  being  dependent  on  the 
temperature  ratio,  will  vary  with  a  variation  of  temperature  loading,  a 
dynamic  similarity  deficiency  may  be  expected  between  cases  of  different 
temperature  loadings.  For  complete  dynamic  similarity  then,  a  temperature 
ratio  must  be  included  in  the  dimensionless  relation,  so  that 

hd  =«  f(Vd  ,  ,  (temperature  ratio).  (3) 

k  k  1 

If  it  is  quite  possible  that  the  nature  of  the  similarity  deficiency 
with  temperature  loading  will  be  different  depending  on  the  direction  of 
heat  flow,  then  two  relations  of  the  above  fora  are  necessary  to  describe 
the  two  processes  -  heating  and  cooling.  In  these  relations  the  fluid 
properties  may  now  be  evaluated  at  any  of  the  temperatures  ,  Ts  } 

or  a  suitable  mean  temperature,  and  the  temperature  ratio  may  be  obtained 
by  any  combination  of  any  two  of  the  temperatures.  In  line  with  some 
of  the  more  successful  previous  investigations,  we  may  choose  the 
arithmetic  mean  temperature  for  the  evaluation  of  all  the  fluid  properties 
(1,  2,  3f  4)  and  a  ratio  of  T  M  and  ?m  for  the  temperature  ratio  (1), 
so  that: 


UK 


IASS Trim 


where  the  subscript  e  refers  to  the  arithmetic  sean  temperature. 

However,  a  correlation  of  the  above  dimensionless  groups 
still  fails  to  give  a  unique  expression  for  the  heat  transfer  process, 
as  is  apparent  by  a  consideration  of  the  following  points! 

a)  Considering  the  fact  that  for  a  particular  fluid  the 
properties  of  the  fluid  do  not  vary  with  temperature  in  the  same  way  at 
different  temperature  levers,  it  say  be  expected  that  the  fashion  the 
group  Tqq  /Tn  enters  in  the  heat  transfer  expression,  i.e.  the  temperature 
loading  factor,  will  be  different  at  different  temperature  levels-  By 
the  same  token  it  will  also  be  dependent  on  the  range  cf  temperatures 
considered.  Thus,  a  correlation  nay  not  hold  for  cases  differing  frssx 
the  case  for  which  the  correlation  was  obtained  with  respect  to  the 
temperature  loading  factor,  and  an  evaluation  of  the  correlation  at  the 
particular  temperature  level  and  range  of  interest  becomes  necessary. 

b)  Apart  from  a  particular  fluid  at  different  temperature 
levels  we  nay  also  consider  a  group  cf  fluids  at  similar  temperature 
levels.  If  the  variation  of  the  transport  properties  with  temperature 
of  the  fluids  considered  are  not  identical,  them  it  becomes  necessary 
to  evaluate  the  temperature  loading  factor  individually  for  each  fluid. 

Since  inclusion  of  the  group  /T  in  the  heat  transfer  expression 

is  to  account  for  the  shift  in  similarity  _»f  the  fluid  property  gradients, 
the  * property  ratios1  may  have  to  be  considered  directly  in  the  correlation, 
in  place  of  the  temperature  ratio,  to  give  a  unique  relation  valid  over 
a  group  of  fluids.  Since  for  a  particular  gas  the  significant  transport 
properties  vary  nearly  identically  with  temperature,  variations  in  all 
the  properties  may  be  lumped  together  by  cons-dering  a  ratio  of  values 
of  any  one  of  them.  Considering  kinematic  viscosity  as  the  property, 
we  obtain  the  mere  general  expression; 

M  =  f(vd,  2^  ,  WM  (5) 

va  .-q- 

1.4  In  light  of  the  above  discussion  ard  the  conditions 
stipulated  in  Section  1.2,  we  cay  now  formalise  the  purposes  behind 
the  present  investigation  as  follows. 

a)  To  investigate  the  nature  of  the  temperature  loading 
factor  for  heating  of  cylinders  in  cross-flow  for  high  temoei ature  loadings. 

b)  To  investigate,  in  terms  of  improved  uniqueness  of  the 
correlation,  whether  it  is  worthwhile  to  replace  the  ratio  **s0  /  V 

/l^j,  when  different  gases  and  gas  mixtures  are  considered  at^high 
temperature  loadings. 

c)  To  investigate  whether  the  slight  variation  in  PTandtl 
nurnb-r,  (due  to  temperature  loading  and  due  to  the  flow  being  of 
different  gases  and  gas  mixtures),  is  significant  enough  to  be 
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considered  in  the  heat  transfer  relation. 

d)  To  obtain  a  general  correlation  for  forced  convective 
heat  transfer  for  the  conditions  stipulated  in  Section  1.2. 

1.5  The  ranges  of  variables  considered  in  the  present 
investigation  were  as  follows: 

a)  Reynolds  number:  The  Rem  (based  on  cylinder  diameter) 
range  considered  was  5-75.  However,  only  data  from  the  range  Rem  ■■  5  to  40 
were  correlated. 

b)  Temperature  loading:  The  maximum  environment  temperature 

(Toq)  considered  was  around  1500°K,  a  limitation  imposed  by  uncertainties 

in  the  evaluation  of  the  fluid  properties  above  that  temperature.  The 

cylinder  temperature  (Ts)  range  considered  was  350°K  to  550°K.  Thus,  the 

temperature  loading  range  achieved  was  Too  s  2  to  4 

T 

l5 

c)  Flow  composition:  The  flows  were  composed  of  N?  and  He 
and  mixtures  of  Ng-He,  N2-CO2  and  He-CC^.  For  each  mixture  at  least  two 
mixture  ratios  were  considered.  However,  onxy  data  from  He  and  N2  flows 
were  considered  for  the  statistical  analysis,  while  the  rest  of  the  data 
were  checked  against  the  result  of  this  analysis.  The  Prandtl  number 
range  amongst  data  for  the  statistical  analysis  was  Pr  >=  0.72  -  0.75. 


2.0  REVIEW  OF  PREVIOUS  INVESTIGATIONS 


2.1  Considerable  data  are  available  for  forced  convective  heat 
transfer  from  heated  cylinders  in  air  cross-flow  (1,2,3,5,6,7*8,9,10). 

On  the  other  hand,  data  are  scarce  for  the  case  of  heat  transfer  to  cooled 
cylinders  (11,  12,  13).  McAdams  (14)  has  collected  and  correlated  results 
from  some  of  these  sources.  However,  due  to  misinterpretation  of  the 
more  important  original  data,  McAdams’  correlation  must  be  considered 
somewhat  approximate  (for  elaboration  on  this  point  see  Ref.  4).  Douglas 
and  Churchill  (4)  recorrelating  some  of  the  data  considered  by  .McAdams 
along  with  some  more  recent  results,  indicated  that  a  unique  correlation 
is  possible  for  cases  of  heating  and  cooling  over  a  wide  range  of  temp¬ 
erature  loadingsby  evaluating  all  the  fluid  properties  in  the  correlation 
at  the  arithmetic  mean  temperature.  However,  considering  the  remaining 
discrepancy  in  results  between  different  investigations  under  dissimilar 
and  even  similar  cases.it  is  to  be  concluded  that  such  a  general  cor¬ 
relation  is  at  best  an  engineering  approximation  and  may  not  be  suitable 
for  accurate  instrument  work. 

Among  the  data  considered  by  Douglas  and  Churchill  (4)  only 
those  of  Hilpert  (2)  and  King  (6)  fall  in  the  low  Reynolds  number  range 
of  our  interest,  i.e.  Rem  <  100.  However,  their  data  are  for  heat 

transfer  from  heated  cylinders.  On  the  other  hand,  all  data  for  heating 
of  cylinders,  i.e.  data  of  Reiher  (12),  Vomhem  (13),  and  Churchill 
and  Briar  (11),  are  for  Reynolds  numbers  above  300.  Still,  the  data  of 
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Churchill  and  Briar  are  of  special  interest  since  they  are  the  only 
ones  obtained  for  any  real  high  temperature  difference  existing  between 
the  cylinder  (100°F)  and  the  gas  stream  (580-1800°F) . 


2.2  Consideration  of  temperature  loading  effect 

Among  the  investigations  reported  by  Douglas  and  Churchill  (4), 
only  those  by  Hilpert  (2)  for  convection  from  cylinders  and  those  by 
Churc’:ill  and  Briar  (11)  for  convection  to  cylinders  consider  explicitly 
the  temperature  loading  effect.  However,  in  the  former  investigations  the 
effect  was  obtained  over  a  narrow  range  of  Reynolds  numbers.  A  consideration 
of  this,  along  with  the  fact  that  Hilpert ’ s  power  law  correlation  is  an 
approximate  representation  of  the  data,  led  Collis  and  Williams  (l)  t-o 
obtain  a  more  accurate  correlation  for  hot-wire  applications.  Further, 
Collis  and  Williams  obtained  their  correlation  under  conditions  of  low 
turbulence  intensity  in  the  free-streanj  the  significance  of  the  effect 
of  free-stream  turbulence  on  heat  transfer  had  not  been  realized  in  previous 
investigations.  The  data  due  to  their  carefully  conducted  experiment  were 
correlated  in  the  following  form: 

^'17=  A+B.Remn  .  (6) 


The  values  of  the  constants  in  the  different  Rejj  ranges  are  shown  in 
Table  I. 


Rem 

A 

’  ' 

B 

n 

0.02-/44 

0.24 

0.56 

0.45 

44-140 

0 

0.43 

0.51 

TABLE  I 


This  may  be  compared  to  Hilpert1 s  correlation: 


Num 


r 


Re, 


ra 


ft) 


0.25 


m 


(7) 


Table  Ilcontains  the  values  of  the  constants  in  the  different  Rejjj  ranges. 


^®m 

c 

m 

1-4 

0.891 

0.330 

4-40 

0.821 

0.385 

40-4000 

0.615 

0.466 

TABLE  II 


The  important  point  to  note  is  that  both  investigations  showed  that 
for  forced  convective  heat  transfer  from  heated  cylinders,  an  increase  in 
temperature  difference  is  associated  with  an  increase  in  hd  .  For  heat 
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transfer  in  the  opposite  direction  ,  i,e.  heat  transfer  to  cooled- 
cylinders,  the  opposite  effect  is  to  oe  expected,  i.e.  hd/km  will  decrease 
with  an  increase  in  temperature  difference.  This  seems  reasonable  since 
an  increase  in  temperature  difference  has  opposite  effects  on  the  velocity 
and  the  temperature  gradients  as  the  heat  transfer  direction  is  reversed. 
Consider  the  case  of  heat  transfer  from  cylinders,  i.e.  Tg  > 

In  this  case  viscosity  of  the  fluid  near  the  wall  is  greater  than  the 
viscosity  in  the  free-stream,  causing  deceleration  of  the  flow  near  the 
wall  which  results  in  a  steeper  velocity  profile  in  the  aerodynamic 
boundary  layer  as  compared  to  the  isothermal  case.  As  the  temperature 
difference  is  increased,  the  viscosity  difference  is  increased,  the  de¬ 
celeration  of  the  flow  near  the  wall  is  greater  and  the  boundary  layer 
profile  is  steeper.  Since  variation  of  thermal  conductivity  with  temperature 
is  similar  to  that  of  viscosity,  a  similar  effect  occurs  for  the  temperature 
profile.  On  the  other  hand,  when  T^,  >  Ts,  the  viscosity  near  the  wall  is 
less  than  in  the  free-stream  -  consequently  the  flow  near  the  wall  acceler¬ 
ates  with  a  flatter  velocity  profi?.e.  As  the  temperature  difference  is 
increased  the  flatness  of  ths  velocity  and  temperature  profiles  will  in¬ 
crease.  From  this  reasoning,  though  no  conclusion  can  be  reached  as  to 
whether  hd/k^  will  either  increase  or  decrease  with  an  increase  of  temp¬ 
erature  difference  for  a  particular  direction  of  he^.t  flow,  it  appears 
clear  however,  that  the  effect  of  a  variation  of  the  temperature  differ¬ 
ence  will  be  opposite  if  the  direction  of  heat  flow  is  reversed. 


However,  Churchill  and  Briar's  (11)  results  from  their 
experiment  with  cooled  cylinders  at  100°?  ir.  heated  nitrogen  flow,  where 
the  Nitrogen  temperature  was  varied  systematically  from  580°?  to  1S00°F, 
show  the  effect  on  hd/l^  with  temperature  loading  to  be  in  the  same 
direction  as  obtained  by  Hilpert  (2)  and  Collis  avid  Williams  (1)  in  their 
experiments  with  heated  cylinders  and  wires  in  ambient  air  flow.  Churchill 
and  Briar's  correlation  is  in  the  following  form: 

NUjjj  =  0.60  Pr^,0’33  ReJ^  3  ^Too^0'12  (8) 


for  a  dG  range  of  300  to  2300. 

/*m 

It  is  this  anomaly  in  the  direction  of  the  temperature  loading 
obtained  by  Churchill  and  3riar  that  focussed  attention  in  the 
investigation  to  be  placed  on  the  determination  of  the  correct 
loading  factor. 


effect 

present 

temperature 


2.3  Consideration  of  the  Reynolds  Number  Dependency 
of  the  Heat  Transfer 


Earlier  investigators,  notably  King  (6),  for  theoretical 
reasons  presupposed  in  their  correlation  the  linear  dependency  of 
Nusselt  number  on  (Rem)^*5$.  Kramer  (5),  analysing  results  from  later 
experiments,  obtained  the  same  1'orm  of  correlation  as  King  but  with 
modified  constants: 

Nu  *  0.42  Pr°'2  -4-  0.57  Pr°*33  Ee°*5 


(9) 
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The  exact  mnfe  of  Re  and  the  reference  temperature  for  procerty  value 
evaluation  were  not  specified.  Cn  the  other  hand,  although  nilpert’s  (2) 
pov;er  lav;  correlation  (Equation  7)  is  somewhat  approximate  and  i(nores 
the  finite  heat  transfer  at  zero  Ren,  it  nevertheless  indicates  that  the 
dependency  of  on  Re,.  is  itself  a  function  of  Re^.  He  noted  disconti¬ 
nuities  in  the  slor-e  of" ‘the  heat,  transfer  curve  at  Re„  =  4  and  Re_  =  40. 

Later,  more  careful  experiments  with  hot-wires  (151  tended 
to  show  a  systematic  deviation  of  observed  Mum  from  the  Num  obtained  from 
King’s  equation  with  increasing  Hem. 

Collis  and  Williams  (l),  provoked  by  this  deviation  from  King’s 
law  and  also  by  the  suspected  discontinuity  at  Ren  =  40  conducted  their 
experiment  very  carefully  and  obtained  a  Re.,  dependency  of  Nu_  different  from 
that  of  King’s  law  and  also  found  a  discontinuity  at  Een  -  44.  This  dis¬ 
continuity  was  confirmed  by  an  aerodynamic  investigation  of  the  flowrfield 
around  the  wire  by  observing  the  onset  of  eddy  shedding  at  Rera  -44  (l). 

Now,  considering  Churchill  and  Briar’s  (11)  analysis  of  their 
data  it  is  to  be  noted  that  they  presupposed  the  linear  relation  of  Nu-,  to 
Ren^'50.  Further,  the  correlation  was  obtained  by  assuming  no  heat 
transfer  at  zero  Eem,  i.e.  they  presupposed  the  form  of  correlation  as: 

0.33  0,5  -  n  ,  , 

Hu,,  =  A  Pr  Rem  ( TuA  (10) 

lT s  / 

and  obtained  the  values  of  the  constants  A  and  n.  Apart  from  the 
assumption  as  to  the  Rera  dependency  of  Num,  this  form  of  correlation  may 
be  justifiable  at  higher  Ben  where  the  heat  f ransx'er  at  zero  Re,^  is  a 
small  percentage  of  the  total  heat  transfer,  but  an  extrapolation  of  the 
correlation  to  low  Be^  may  not  be  accurate  enough.  It  must  be  noted  that 
their  correlation  was  obtained  for  dG/um  beyond  300. 

Thus,  in  the  present  investigation  at  low  Fe_j,  the  determination 
of  the  exact  Re^  dependency  of  Num  is  considered  as  part  of  the  analysis. 

2.4  Consideration  of  Brandt!  Humber  Effect 

Kramer’s  (5)  analysis  explicitly  considers  the  effect  of 
Prandtl  number  variation  on  heat  transfer.  Data  were  obtained  for  air 
(Pr  -  .71),  water  (Pr  *  7.1),  white  oil  (Pr  -  29.6),  Transformer  Oil  I 
(?r  =*  35.)  and  Transformer  Oil  II  (Pr  »  525.).  From,  this  correlation 
(see  equation  9)  it  appears  'hat  if  interest  is  limited  to  heat  transfer 
with  gases  only,  the  effect  of  the  slight  differences  in  Pr  for  different 
species  v?ill  be  very  small  and  may  not  be  detectable  in  an  experiment,  due 
to  normal  data  scatter.  By  the  same  token,  the  effect  of  the  small 
variation  of  Pr  with  temperature  for  a  particular  species  of  fas  will  also 

f{i  be  negligible. 

Ill 

Investigations  with  different  gases  have  not  been  carried 
out,  all  investigators  having  used  air  with  the  exception  of  Churchill 
and  Briar  (11)  who  used  Nitrogen.  As  such,  it  is  not  possible  to  determine 
the  uniqueness  of  any  temperature  loading  factor  for  different  gas  species  - 
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a  point  to  be  investigated  in  the  present  analysis. 


3.0  DESCRIPTION  OF  EQUIPMENT.  EXPERIMENTAL  SET-UP 

AND  MEASUREMENT  PROCEDURE 

3-1  Heat  transfer  measurements  were  obtained  by  traversing 
internally  cooled,  electrically  heated,  constant  temperature,  cylindrical 
film  sensors  across  the  potential  core  of  a  plasma  jet.  The  power  input 
to  the  plasma  jet  and  the  plasma  gas  flow  rate  was  maintained  such  that 
not  only  would  ionization  be  negligible  at  the  jet  exit  plane  but 
also  recombination  of  the  dissociated  species  would  be  complete  in  the 
plenum  chamber.  Binary .mixture  flows  were  obtained  by  mixing  the  second 
gas  specieswith  the  primary  ionization  gas  in  the  plenum  chamber.  Temper¬ 
ature  and  velocity  measurements  were  carried  out  by  means  of  a  Pt  -  Pt 
10%  Rh  thermocouple  and  a  pitot  probe  respectively,  both  mounted  on  the 
same  traverse  as  the  cooled-film  sensor.  All  measurements  as  functions 
of  the  traverse  distance  along  the  diameter  of  the  jet  were  recorded  on 
an  X-Y  plotter.  Knowing  the  relative  position  of  each  instrument  with 
respect  to  the  jet  axis,  point  measurements  of  pitot  pressure,  thermo¬ 
couple  output  and  the  film  sensor  bridge  voltage  were  obtained  from  the 
recorded  prefiles. 

3.2  Description  of  Special  Equipment 

The  sensors,  probes,  electrical  circuitry  and  the  sensor 
ccoling  water  supply  system  are  manufactured  by  the  Thermo  Systems  Inc . , 
Minneapolis,  Minnesota,  as  standard  equipment.  The  description  of  the 
equipment  that  follows  has  been  obtained  from  private  communication  k : th 
Dr.  L.  Fingerson  of  the  Company  and  from  their  various  bulletins,  technical 
notes  and  instruction  manuals,  listed  as  Ref  16. 

3-2.1  Cooled  Film  Sensor  and  Sensor  Holder 

A  photograph  of  the  sensor  with  explanatory  sketches 
is  shown  in  Fig.  2.  The  loop  type  sensor  is  made  of  vycor  tube.  .015 
cm.  0. D.  and  .010  cm-  I.D.  A  thin  film  of  platinum,  1GC0  -  2CC0  8 
thick,  on  one  limb  cf  the  loop  acts  as  the  resistance  element  and  is 
isolated  b7  heavy  gold  plating  (1^-1  mil  thick)  on  the  rest  of  the  loop. 
The  length  of  the  sensing  element,  in  the  sensors  used,  i3  .103  cm  and 
its  resistance  around  3  ohms.  The  sensing  element  is  coated  with  a  film 
(5000  8  thick)  cf  silicon  monoxide  to  reduce  catalytic  effects  when  the 
sensor  is  exposed  to  flame  environment. 

The  sensor  is  soldered  onto  the  tip  of  the  sensor 
holder,  which  serves  as  &n  electrical  lead  and  also  as  the  sensor  cooling 
water  connection  and  contains  a  built-in  filter  on  the  water  inlet  side. 

The  fittings  for  the  electrical  connection  to  the  power  supply  system 
and  the  fittings  for  the  coolant  line  connection  to  the  water  supply 
system  are  mounted  on  a  fitting  plate,  as  shown  in  the  photograph  in 
Fig.  3A. 
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3.2.2  Probe 

A  photograph  and  a  diagram  of  the  probe  with  the 
sensor  holder  assembly  insert oq  are  shown  in  Fig.  38  and  in  Fig.  4 
respectively.  A  close-up  oi  the  probe  tip,  showing  the  sensor  in 
operating  position,  is  shown  in  Fig.  5.  The  probe  is  made  up  of  a 
series  of  cascading  cooling  jackets,  the  cooling  jackets  in  turn  made 
up  of  a  tightly  fitting  bundle  of  stainless  steel  tubes,  to  r.r»ovide 
uniform  and  near  optimal  cooling.  The  probe  is  provided  witi.  . 
traversing  sleeve  and  a  fitting  plate  containing  the  fittings  for  the 
jacket  cooling  water  connections.  When  the  sensor  holder  is  in  place, 
its  fitting  plate  is  attached  to  that  of  the  probe  by  means  of  screws 
in  slots,  allowing  for  a  certain  amount  of  play,  so  that  the  holder  can 
be  pulled  back  and  the  sensor  withdrawn  inside  the  probe  when  it  is  not 
in  operation. 


3.2.3  Sensor  Cooling  Water  Supply  System 

The  sensor  water  supply  system  supplies  coolant  to 
the  sensor  at  constant  pressure  and. hence,  at  a  constant  mass  flow  rate. 

It  consists  of  a  stainless  steel  pressure  tank  of  1  litre  capacity  which 
is  filled  with  distilled  water  and  pressurized  through  a  regulator  from 
a  nitrogen  bottle. 

3.2.4  Probe  Jacket  Water  Supply  System 

Ordinary  tap  water,  filtered  through  a  25  micron 
filter,  is  used  for  probe  cooling.  The  water  pressure  i3  maintained 
arcund  80  p.s.i.g. 

3.2.5  Bridge  and  Control  Circuit 

A  simplified,  functional  schematic  diagram  of  the 
circuit  is  shown  in  Fig.  6.  The  platinum  film  on  the  sensor  is  maintained 
at  constant  resistance  and  hence  at  a  constant  temperature  by  placing 
it  across  one  arm  of  a  bridge.  The  desired  resistance  of  the  censor 
(calculated  from  the  desired  operating  temperature  and  the  cold  resist¬ 
ance  of  the  sensor)  is  obtained  by  adjusting  the  resistance  deck  Rl. 

A  choice  of  two  bridges  is  possible,  one  for  low  energy  dissipation 
(.005  to  approximately  0. 5  watts),  as  encountered  in  usual  hot-wire  or 
hot-film  operation,  and  the  other  for  higher  energy  dissipation  (up  to 
40  watts)  for  cooled-film  operation. 

Once  a  voltage  is  applied  across  the  bridge  and 
the  bridge  is  balanced,  the  film  will  be  maintained  at  the  desired 
temperature.  Any  change  in  environmental  conditions  (changes  in  temper¬ 
ature,  velocity  or  composition)  causing  a  variation  in  the  rate  of  heat 
transfer  to  or  from  the  film,  will  tend  to  vary  tne  film  temperature  and 
hence  the  film  .-ccistance,  causing  an  imbalance  between  the  two  sides  of 
the  bridge.  This  imbalance  is  amplified  by  the  bridge  amplifier,  which 
signals  the  power  amplifier  to  increase  or  decrease  the  power  to  the 
sensor  according  to  the  direction  of  voltage  change  in  the  bridge.  This 
fsed-back  loop  has  a  frequency  response  of  500  KHz#  (nominal). 
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The  system  allows  measurement  of  both  bridge  voltage 
and  actual  electrical  power  supplied  to  the  film.  The  latter  is  obtained 
by  means  of  a  power  computer  which  takes  values  of  the  bridge  voltage  and 
sensor  resistance  and  uses  a  squaring  circuit,  amplifier  and  a  voltage 
divider.  The  panol  meter  may  be  switched  either  to  read  the  mean  bridge 
voltage  or  the  mean  energy  dissipation  in  watts,  when  the  sensor  is  ex¬ 
posed  to  a  steady  state  environment.  A  photograph  of  the  panel  along  with 
that  of  the  sensor  water  supply  system  is  shown  in  Fig.  7- 

A  detailed  description  of  the  circuit  may  be  found 

in  Ref.  16A. 


3*2.6  Cooled-rilm  Sensor  Operating  Principle 

Referring  to  Fig.  1  we  note  that  for  usual  constant 
temperature,  cooled-film  application,  i.e.  T^  >  T  ,  a  heat  balance  of  the 
system  gives : 

^supplied  ”  qcoolant  ~  %sv  - 


If  the  cooling  water  inlet  temperature  and  flow  rate 
(i.e.  the  coolant  pressure)  is  kept  constant,  QqoclANT  re2a^n  constant. 

A  measure  of  ^SUPPLIED  then  will  be  a  measure  of  wjjjy,  i.e.,  the  heat 
transferred  from  the  environment  to  the  sensor. 

For  adequate  operation  it  is  then  mandatory  that 
QcooLANT  >  OSKV  (MAXIMUM)*  LI  heat  transferred  to  the  coolant  is  less 
than  that  absorbed  from  the  environment,  not  only  does  measurement  bet* see 
impossible  but  the  sensor  temperature  will  increase  beyond  the. set  value 
and  the  sensor  may  burn  out.  A  proper  QcOGLAKT  s*/  he  set  by  exposing 
the  sensor  to  a  condition  of  maximum  heat  transfer  and  then  by  adjusting 
either  the  coolant  pressure  or  the  sensor  temperature  or  both  to  obtain 
a  Qcoolant  greater  than  Qsjjy  (MAX)-  Since  a  prior  determination  of  (MAX) 
is  not  always  possible  (viz.  as  in  the  case  of  application  in  the  hyper¬ 
sonic  wakes  of  projectiles  launched  in  free-flight  ranges),  a  rough  estimate 
of  the  expected  maximum  heat  transfer  may  be  made,  and  from  the  calibration 
chart  in  Fig.  8,  jthe  sensor  operating  temperature,  and  coolant  pressure  may 
be  determined.  As  shown  in  Fig.  8,  at  a  particular  coolant  pressure  the 
maximum  sensor  temperature  is  limited  because  of  coolant  boiling  inside 
the  sensor  resulting  in  unstable  heat  transfer  to  the  coolant. 


The  sensor  is  set  to  operate  at  the  desired  temp¬ 
erature  by  setting  the  proper  resistance  in  the  resistance  deck.  The 
operating  resistance  is  obtained  by  measuring  the  sensor  '..old"  resist¬ 
ance  and  using  the  relation: 


R 

s 
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where  Rg  is  the  sensor  resistance  at  the  selected  operating  temperature, 

T  . and  ?•<.  is  the  "cold"  resistance  at  temperature  ?c  (usually  room  temper¬ 
ature)  .  The  temperature  coefficient  of  resistance/  of.* is  obtained  by 
a  calibration  of  a  representative  number  of  sensors  or  fros  the  manufacturer* 
quoted  figures.  A  typical  calibration  curve  for  of  obtained  for  three 
sensors  free  the  batch  used  in  the  present  investigation  is  shown  in 
?ig.  9- 


The  value  of  QcqOLAHt  measured  by  maintaining  the 
sensor  at  the  operating  temperature  and  coolant  pressure  and  exposed  to 
room  environment.  Since  conservatively  calculated  heat  losses  from  the 
sensor  due  tc  radiation,  free  convection  and  conduction  together  add  up 
to  be  two  orders  of  magnitude  less  than  the  heat  transferred  to  the  coolant; 
ail  the  heat  transferred  under  room  condition  cay  be  considered  to  be 
transferred  to  the  coolant,  i-e. 


n  *  Q 

SUPPLIED  (at  room  temp. )  COOLANT 


QcODLAKT  say  be  measured  directly  in  electrical  units 
by  using  the  power  computer "or  by  measuring  the  bridge  voltage  and  using 
the  following  relation  (refer  to  Fig.  6): 


n 


where 


5 

3 


Zt  a 
t  s 


(•O  .  P  .  3 

's'S  *'d  ^  Tr«ra/ 


‘3  *  “CCK?' 
bridge  voltage 
sensor  resistance 

resistance  of  the  top  leg  of  the  bridge 
probe  and  cable  resistance. 


For  the  measure  of  Qgnspr rrn  under  operation  the 
same  two  methods  may  be  used-  However,  since ^tHere  was  some  doubt  as 
to  the  calibration  of  the  power  computer,  in  the  present  investigation 
it  was  the  bridge  voltage  that  was  recorded. 


3-3  Experimental  Set-Up 


A  schematic  and  a  photograph  of  the  experimental  set-up  aTe 
shown  in  Figs.  10  and  11  respectively.  The  hot  jet  was  obtained  by  means 
of  a  bench  counted  plassadyne  plasma  torch  fitted  with  a  nozzle,  the 
nozzle  exit  diameter  being  $  in.  The  gas  flow  rates  were  metered  by  means 
of  calibrated  flow  meters  and  pressure  gauges  mounted  on  the  control 
panel  of  the  plasma  torch.  The  measuring  instruments ,  i.e.  the  cooled 
fils  sensor  probe,  the  thermocouple  and  the  pitot,  were  mo  untea  on  a 


UNCLASSIFIED 

12 


traversing  table.  A  close-up  photograph  of  the  instruments  and  the  plasma 
torch  is  shown  in  Fig.  12.  The  traverse  distances  of  the  table  along  and 
across  the  .jet  were  monitored  by  a  knife-edge  slider  fixed  to  the  table 
and  sliding  over  a  calibrated  linearly  variable  potentiometer.  The  cross- 
traverse  potentiometer  output  was  connected  to  the  fX’  in  put  of  a  Hosley 
X-Y  plotter,  whereas  the  axial  traverse  potentiometer  reading  was  displayed 
on  a  John  Fluke  differential  vacuum  tube  voltmeter. 

A  ?t  -  Pt  10?  Rh  thermocouple  with  a  bead  diameter  of  .127  cm 
was  used  to  measure  temperature.  It  was  thought  to  be  more  judicious  to 
use  an  unshielded  thermocouple  and  to  correct  for  radiation  and  conduction 
loss  rather  than  to  use  a  shielded  thermocouple  with  uncertain  radiation 
effect.  Further,  the  minimum  size  limitation  and  the  vulnerability  of  the 
shielding  to  high  temperature  with  consequent  geometry  change  precluded  the 
use  of  shielded  thermocouples. 

For  velocity  measurements  a  carbon-tipped  pitot  probe  with  an 
orifice  sice  of  .103  cm  was  used.  The  pitct  pressure  was  measured  by  means 
of  a  calibrated  HXS  Baraton  type  77  pressure  aster  using  a  type  774 
pressure  head. 

The  thermocouple  mv  output,  the  pressure  meter  output  and 
the  bridge  voltage  of  the  film  sensor  were  connected  to  the  TI*  input  of 
the  X-I  plotter  through  a  control  box,  so  that  during  a  run  as  each 
instrument  would  traverse  across  the  jet  their  outputs  could  be  successively 
switched  onto  the  *TT  input  of  the  recorder  for  recording.  The  recording 
station  arrangement  is  shewn  in  Fig.  13. 

The  cross  traverse  distance  was  recorded  at  double  scale  for 
better  positional  resolution,  whereas  the  outputs  of  each  instrument  were 
recorded  on  a  scale  to  give  the  maximum  deflection  on  the  recorder.  A 
typical  recording  of  a  run  is  shown  in  Fig.  It. 

3-4  Measurement  Procedure 

3.4.1  Instrument  Alignment 

Before  each  set  of  runs  (usually  consisting  of  4  to  5 
runs)  the  sensing  elements  of  the  measuring  instruments  were  aligned  to 
be  in  ths  same  horizontal  and  vertical  planes,  the  horizontal  plane  being 
that  of  the  jet  axis  and  the  vertical  plane  that  of  the  cross-traverse, 
i.e.  normal  to  that  of  the  jet  *x±s.  The  eooied-fiis  sensor  was  aligned 
to  be  accurately  normal  to  the  jet  axis.  The  alignment  of  the  instruments 
was  achieved  by  a  microsccpe  replacing  the  plasma- jet  on  its  bracket 
such  that  the  axis  of  the  jet  would  coincide  with  that  of  the  microscope, 
as  shewn  in  Fig.  15,  and  by  a  theodolite  placed  at  a  slight  angle  to  the 
vertical  plane  of  the  instruments.  Alignment  was  ensured  by  traversing 
the  instruments  successively  to  the  microscope  axis  and  adjusting  their 
positions  such  that  they  coincided  with  the  cross-hairs  of  both  the  micro¬ 
scope  and  the  theodolite.  After  alignment  the  positions  of  the  instruments 
with  respect  to  the  jet  axis  and  each  other  were  noted  by  traversing 
them  to  the  microscope  axis  and  recording  the  cross-traverse  potentio¬ 
meter  output  on  the  X-I  plotter. 
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During  alignment  the  physical  condition  of  the 
sensing  elements  as  to  deposition  of  contaminant  on  the  cooled  film, 
flaking  of  the  film  and  thermocouple  geometry-  alteration  was  checked. 

3-4.2  Plasma  Jet  Conditions 

Power  Setting 

The  plasma  jet  was  run  at  very  low  power  inputs  so 
that  the  temperature  at  the  no2zle  exit  never  exceeded  2000°K.  This  was 
to  ensure  that  no  part  of  the  gas  flow  was  ionized.  Further,  the  gas  flow 
rates  were  such  that-  the  residence  time  in  the  plenum  chamber  was  sufficient 
for  the  recombination  of  the  products  of  dissociation.  (The  latter  con¬ 
dition  was  verified  by  the  addition  of  CO^  in  the  plenum  chamber  and  noting 
the  absence  of  any  CO  flame  at  the  jet  exit.) 

Jet  Stability 

The  stability  of  the  jet  after  a  5  min.  warm-up  period 
was  checked  by  placing  the  cooled-film  sensor  in  the  jet  and  recording  the 
bridge  voltage  on  a  time  scale.  Since  the  film-  will  sense  any  change  in 
both  the  velocity  and  the  temperature  of  the  jet,  this  ensured  a  check 
cm  any  drift  of  both  properties.  Over  a  recorded  period  cf  15  minutes 
(average  time  for  a  run)  no  drift  was  noted.  However,  free  time  to  time 
the  .let  would  tend  to  become  unstable  due  to  electrode  erosion.  Headings 
obtained  during  these  periods  were  either  discarded  or  retaken. 

Jet  Turbulence 


The  turbulence  level  in  the  core  of  the  jet  was 
estimated  by  monitoring  the  cooled-film  bridge  voltage  by  an  oscilloscope. 

The  maximum  cooled-film  frequency  response  is  around  50  5.C.S.  Since  the 
film  will  sense  not  only  the  velocity  fluctuations  but  the  temperature 
and  concentration  fluctuations  as  well,  decoupling  cf  the  signal  to  obtain 
quantitative  measurement  of  the  intensity  of  each  parameter  becomes  cos- 
plicated.  Free-  the  oscilloscope  tracing,  a  tracing  of  the  superimposed 
effects  of  the  fluctuation  of  all  the  three  parameters  in  the  core  of  the 
jet,  the  maximum  relative  intensity  was  found  to  be  less  than  10%  ^  fil*  01). 
If  it  is  considered  that  the  film  is  sensing  only  velocity  “v 

fluctuations,  the  most  acute  situation,  then  the  error  in  heat  transfer 
measurement  will  be  less  than  Z%.  (See  Ref.  1?  for  an  elaboration  cm 
the  effect  of  large  turbulence  fluctuaticns  on  heat  transfer. )  The 
error  due  to  temperature  fluctuation  will  arise  because  of  consideration 
of  the  fluid  properties  evaluated  at  the  average  temperature  rather  than 
a  consideration  of  the  average  of  the  fluid  properties  -  the  difference 
between  then  being  due  to  the  non-linear  variation  of  fluid  properties 
with  temperature.  Similarly  the  error  because  of  concentration  fluc¬ 
tuations  viil  be  due  to  the  difference  between  the  evaluated  properties 
at  the  mean  concentration  rather  than  the  average  of  the  properties  -  a 
difference  because  cf  non-linear  variation  of  fluid  properties  with 
percentage  concentration.  However,  these  effects  will  be  sail  compared 
to  these  duet©  large  velocity  fluctuation,  i.e.  due  to  the  lateral  tur¬ 
bulence-velocity  c  amponent  perpendicular  to  the  cylinder. 


UNCLASSIFIED 

14 


3-4.3  Run  Conditions 

During  stable  operation  of  the  jet,  the  pitot,  the 
thermocouple  and  the  cooled-fiim  sensor  were  traversed  across  the  jet,  the 
last  traversing  thrice  for  three  different  sensor  resistances  and  hence 
three  different  temperature  settings.  Due  to  the  variatio/i  of  the  sensor 
“cold"  resistance  from  one  sensor  to  another  and  sometimes  from  one  run  to 
another  it  was  not  always  possible  to  set  the  sane  three  sensor  temperatures 
for  all  runs.  The  distances  of  the  section  of  traverse  from  the  jet  exit 
were  such  that  the  major  portion  of  the  traverse  would  be  through  the 
potential  core  of  the  jet.  A  typical  recording  of  a  single  run  is  shown  in 
Fig.  14.  To  obtain  measurements  over  a  range  of  ReE  and  temperature  loading 
the  flow  rate  of  the  cas  and  the  plasma  power  input  were  varied.  For  each  gas 
species  or  species  mixture,  usually  three  power  settings  and  for  each 
power  setting  three  flew  rates  were  employed.  With  a  particular  power 
setting  and  flow  rate,  usually  three  or  four  traverses  at  different  distances 
from  the  jet  exit  plane  were  made.  For  mixture  runs,  three  different 
mixtures  -  S2  -He,  !<2  -CO2,  He  — CO2  were  considered.  During  these  runs, 
usually  the  major  species  interms  of  mol  fraction  was  used  as  the  arc-gas, 
the  second  speciesbeing  mixed  with  this  in  the  plenum  chamber.  However, 

CO2  was  never  used  as  the  arc-gas.  For  each  mixture,  two  mixture  ratios 
were  considered.  Due  to  the  volume-flow  rate  limitation  of  the  plasma-jet. 
judicious  settings  for  the  mixture  ratio  were  not  always  possible. 


4.0  DATA  ANALYSIS  PROCEDURE.  RESULTS  AND  DISCUSSION 

4.1  The  values  of  thermocouple  eaf,  pitot  pressure  and  the  cooled  - 
film  bridge  voltage  at  any  particular  point  in  the  traverse  were  obtained 

by  taking  values  equidistant  from  the  jet  axis  from  their  respective  profiles, 
the  position  of  the  jet-axis  with  respect  to  each  traverse  having  been  pre¬ 
recorded  as  mentioned  in  Section  3.4.1.  Such  point  values  were  obtained 
free  only  the  potential  core  of  the  jet. 

4.2  The  unconnected  temperature  was  obtained  from  the  millivolt 
reading  of  the  thermocouple  using  the  calibration  tables  of  Omega  Eng.  Inc. 
This  temperature  was  corrected  for  the  error  due  to  radiation  from  a  heat- 
balance  of  the  thermocouple  bead.  (See  any  standard  text  for  elaboration). 
Free  an  independent  approximate  calculation  the  conduction  error  was  found 
to  be  negligible.  For  the  heat  balance,  the  thermocouple  bead  was  considered 
to  be  a  sphere  and  the  forced  convective  heat  transfer  coefficient  was 
obtained  from  Kramer* s  ( 5)  correlation  for  spneres,  and  the  emissivity  factor 
for  the  pi  at. inns  bead  as  a  function  of  temperature  was  obtained  from  Ref.  IS. 
Since,  evaluation  of  the  forcea  convection  coefficient  demands  a  determination 
of  flow  velocity  and  the  fluid  transp'rt  properties  at  the  true  temperature 
and  since  tne  velocity  interpretation  from  the  pitot  pressure  is  again 
dependent  on  the  knowledge  of  fluid  density  at  the  true  temperature,  an 
iteration  process  was  used  to  solve  for  the  true  temperature. 

The  transport  property  values  of  the  gas  spdeies  and  species 
mixtures  for  the  above-mentioned  temperature  correction  process  and  for 
the  subsequent  data  analysis  were  calculated  from  the  expressions  collected 
in  Ref.  19.  Details  of  the  expressions  used  may  be  found  in  Appendix  A. 


UNCLASSIFIED 

15 


From  the  coojed-film  bridge  voltage  outside  the  jet  and  from 
the  bridge  voltage  when  fils  was  at  the  point  of  measurement  in  the 
jet,  the  heat  transferred  from  the  environment  was  obtained  according  to 
the  method  discussed  in  section  3-2.0.  The  heat  transfer  coefficient  was 
then  determined  frear 


is/ 

4-3  As  a  preliminary  step  to  data  correlation  the  dimensionless 
heat  transfer  coefficient  hd/k_j  was  plotted  3s  a  function  of  He_  to  determine 
the  general  trends  and  the  data  scatter.  The  following  points  were  noted. 
(Fig.  26  ray  be  consulted  for  this  discussion.) 

(i)  Because  of  the  narrow  range  of  volume-flow  rates  possible 
in  the  plasma- jet,  the  Re^,  range  obtained  for  a  particular  species  or  species 
mixture  by  a  variation  of  the  plasma-gas  flow  rate  was  rather  limited. 
However,  because  of  the  widely  differing  kinematic  viscosities  of  the  species 
and  species  mixtures  considered, the  overall  Rs^  range  obtained  extended 
from  5  to  75- 

(ii)  If Ujj,  data  for  Re-  more  than  55  were  above  the  trend  line 
of  the  rest  cf  the  points  (i.e.  those  below  Re-  «  40}when  extrapolated  to 
that  range.  Ho  data  were  obtained  for  R©q  between  40  and  55-  From  previous 
investigations  (1,  2)  with  heated  cylinders,  a  change  in  slope  in  the  heat 
transfer  curve  was  noted  around  Re-  «  40  -  44  due  to  the  onset  of  eddy 
shedding.  Sowever  there  is  no  reason  to  expect  the  onset  of  eddy  shedding 
at  exactly  the  same  Re^  for  the  case  of  cooling  of  cylinders  since  the 
boundary  layer  characteristics  in  the  two  cases  are  different,  as 
discussed  in  section  2.1.  No  aerodynamic  investigation  specifically  aimed 
at  the  determination  of  the  onset  of  eddy  shedding  was  carried  out  in  the 
present  investigation  and,  hence,  any  conclusion  as  to  the  exact  ?Jt&  of 
discontinuity  in  the  heat  transfer  curve  and  as  to  the  reason  for  the 
higher  Hu_  for  Re-  greater  than  55  will  be  deferred  until  further  ex¬ 
perimentation.  For  the  d resent  analysis  the  range  of  Ee_  considered  was 
*5  to  4C-- 

(iii)  The  scatter  of  data  ires  the  gas  mixture  runs  was 
noted  to  be  slightly  more  than  that-  f ice  the  pare  gas  runs.  Possibly 
this  is  due  to  imperfect  mixing  of  the  gases  in  the  plenum  chamber  of 
the  jet.  For  this  reason  and  also  because  of  less  certainty  in  the 
property  value  evaluation  of  mixtures  at  high  temperatures,  the  subsequent 
statistical  analysis  was  confined  to  data  from  pare  He  and  flews  only, 

(iv)  Due  to  the  narrow  ran--?  of  Re,,  obtained  for  a  particular 
gas  species  or  species  mixture  and  due  to  the  data  scatter,  the  expected 
slight  influence  of  ths  very  small  variation  of  Prandtl  number  on  heat 
transfer  could  not  be  properly  detected.  In  any  case,  since  the 
statistical  analysis  was  to  be  confined  to  data  from  He  and  Kj  flows,  the 
Prandtl  numbers  being  very  nearly  identical  for  the  two  gases,  con¬ 
sideration  of  Frandti  number  as  a  variable  was  not  considered  important. 
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4.4  In  the  light  of  the  above  discussion,  it  was  decided  to 
correlate  the  data  from  pure  He  and  N2  flows,  between  Rem  »  5  and  Rem  =  40, 
by  a  relation  of  the  following  form: 


Num  (T ~a_y 


=  A1 4-  B1 


Re 


mJ 


m 


the  values  or  the  constants  n~,  A^,  B  and  m  to  be  determined  by  a  least 
square  fit.  The  values  of  the  constants  were  evaluated  separately  for 
N2  and  He  flows  to  determine  whether  the  temperature  loading  effect  was 
different  for  the  two  cases.  Further,  the  values  were  evaluated  individually 
for  the  two  cases  of  temper/  ture  loading ,  once  when  the  loading  was  applied 
by  a  variation  of  Tg  while  T  <*,  was  constant  and  again  when  T  ^  was  varied 
while  Ts  was  kept  constant.  This  was  done  to  determine  whether  a  weighted 
mean  temperature  rather  than  the  arithmetic  mean  temperature  was  more 
suitable  for  the  evaluation  of  fluid  property,  however,  since  for  the 
individual  cases  the  Rejjj  ranges  were  limited,  the  value  of  the  constant  m-*- 
w»s  evaluated  by  considering  all  data  together,  and  using  this  value  of  nrA 
the  values  of  r.  ,  A^,  were  calculated  for  each  case.  The  result  of 
this  analysis  is  shown  in  the  table  of  Fig.  16.  The  n1  values  were  found 
to  be  somewhat  different  for  the  two  cases,  though  for  a  particular  gas 
the  values  were  close  for  the  two  types  cf  temperature  loading. 

4.5  The  next  step  was  to  rep?  ice  t,he  temperature  ratio  by  a 
kinematic  viscosity  ratio  in  the  analysis,  so  that  the  correlation  form 
now  became: 

Num  (  —  T'  =  A  ■+■  B .  Re  m 


The  values  of  the  constants  evaluated  fo.-  each  case  and  when  all  He  and 
N2  data  were  considered  together  are  shown  in  the  table  cf  Fig.  17.  Com¬ 
pared  to  values,  the  values  of  n  were  found  to  be  close  for  the  different 
cases.  It  remained  now  to  consider  whether  a  unique  relation  wav  possible* 
The  fits  of  the  data  from  each  case  were  checked  with  the  overall  relation 


whose  constants  were  evaluated  by  considering  all  He  and  data  togethe- 
and  as  shown  in  the  table  of  Fig.  17.  The  worst  increase  in  r.m.s.  error 
was  less  than  2.5$. 


The  results  of  this  analysisare  shown  graphically  in  Figs. 
18  to  26.  Figs  18A  and  19A  are  Num  vs.  Rem  plots  c  0  the  data  from  N2  and 
He  flows  respectively  for  the  temperature  loading  obtained  by  a  variation 
of  T^  while  Ts  was  kept  constant. 

In  Figs.  18B  and  198  are  shown  the  same  data  but  with  the 
temperature  loading  factor  included.  Fig.  20A  and  21A, and  20B  and  21B 
are  similar  plots  but  for  data  obtained  when  the  temperature  loading  was 
applied  by  keeping  T  ^  constant  while  Tg  was  varied. 
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The  next,  step  was  to  consider  the  temperature  loading  effect 
on  data  from  the  mixture  flows.  Considering  the  lower  accuracy  in  the 
evaluation  of  mixture  transport  properties,  it  was  felt  that  the  signif¬ 
icance  of  the  T. L.F.  obtained  for  the  mixture  flows  individually  would 
be  ambiguous  and  in  any  case  not  more  significant  than  the  value  already 
obtained  for  He  and  N  flows.  As  such  the  data  from  the  mixture  runs 
were  checked  for  the  elimination  of  the  temperature  loading  effect 
by  using  the  T.L.F.  obtained  for  He  and  N2  flows.  Figs.  22A,  23A  and 
24A  show  Niijjj  data  from  He-N2,  N2-CO2  and  He~CC>2  mixture  flows  respectively 
(two  mixture  ratios  considered  for  each  mixture)  plotted  as  a  function  of 
Rem.  Figs  22B,  23U  and  24B  are  replots  of  the  corresponding  data  but 
considering  the  T.L.F.  A  considerable  improvement  in  data  scatter  may  be 
noted,  though  no  claim  to  the  complete  elimination  of  the  temperature  loading 
effect  can  be  made. 

All  data  corrected  for  the  temperature  loading  effect  plotted 
together  are  shown  in  Fig.  25-  It  is  evident,  tnat  except  for  CO2-N2  data 
which  were  for  Rem  greater  than  55,  data  below  Rem  -  40  overlap  and  follow 
the  same  heat  transfer  curve  ar.d,  as  such,  show  that  along  with  a  unique 
T.L.F.  a  unique  correlation  is  also  possible.  Finally,  the  obtained  Re,,, 
dependancy  of  the  Num  was  checked  for  the  different  species  mixture  flows 
by  plotting  the  data  oil  Num(  vs.  Rem-^5  grid,  as  shown  in  Fig. 

26.  The  r.m.s.  deviation  of  Num  from  the  correlation  line  was  found  to 
be  .0924  when  all  data  below  Rera  =  40  were  considered  together. 

4.6  Comparing  the  correlation  obtained  in  the  present  in¬ 
vestigation,  i.e.: 


Nuj,,  =  (0,2068  4- .4966  Rem0’45)  (^/^T*15  , 


with  previous  correlations  (see  section  2.0),  the  following  points  may 
be  noted. 


(i)  It  is  found  that  for  the  case  of  heat  -transfer  to 
cooled  cylinders,  as  the  temperature  loading  is  increased,  the  hd/kjjj  is 
decreased.  This  is  consistent  with  the  previously  observed  (1,  2) 
increase  in  hd/km  with  increasing  temperature  loading  for  the  case  of 
heat  transfer  from  heated  cylinders,  as  has  been  elaborated  in  section  2.2. 
However,  from  Churchill  and  Briar’s  (11)  result  for  similar  direction  of 
heat  flow, as  in  the  present  investigation, the  opposite  effect  is  noted. 

(ii)  In  the  Rem  range  of  5  to  40,  the  Num  was  found  to  vary 
linearly  with  Rera*^5.  This  deviation  from  King’s  law  is  similar  to  that 
observed  by  Collis  and  Williams  (1)  in  their  experiment  with  heated  wires. 
However,  there  is  no  reason  to  expect  the  same  law  to  hold  for  the  two 
cases  of  cooling  and  heating  since,  as  previously  discussed,  the  aero¬ 
dynamic  features  in  the  two  cases  are  different.  A  valid  comparison 
exists  only  with  correlations  obtained  for  similar  direction  of  heat 
flow.  However,  Churchill  and  Briar’s  analysis  presupposed  a  linear 
Num  variation  with  Rem- 5  and  as  such  can  not  be  properly  compared  to  the 
present  result. 
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(iii)  Finally  Fig.  2?  compares  values  of  Nun  calculated  from 
some  of  the  correlations  obtained  for  heat  transfer  from  heated  cylinders 
with  values  calculated  from  the  present  correlation,  the  evaluation  being 
for  an  identical  condition  of  heating  of  cylinder  with  T-  -  400°K,  = 

1200°K  ajic  e  flow  composed  of  N2.  Even  considering  the  wide  discrepancies 
in  values  calculated  from  correlations  obtained  for  cooling  of  cylinders 
only,  the  error  involved  in  using  them  to  predict  heat  transfer  for  heating 
of  cylinders  is  quite  apparent. 


5.0  CONCLUSIONS 

5.1  For  forced  convective  heat  transfer  to  circular  cylinders 
placed  normal  to  a  heated  flow,  a  correlation  was  obtained  statistically 
from  data  measured  in  the  Reynolds  number  range  of  5  to  40,  temperature 
loading  range  of  Tco  /Ts  *  2  to  4  and  for  flows  of  N2  and  He.  Over  the 
entire  range  the  following  correlation  was  determined; 

Num  (  )‘i5  =  .2068+  .4966  Rem'45  y 

/>n. 

where  all  the  fluid  properties  in  the  dimensionless  parameters  are 
evaluated  at  the  arithmetic  mea  temperature. 

5.2  Because  of  increased  data  scatter  and  also  because  of  some 
uncertainty  in  the  property  value  evaluation  of  mixtures  at  elevated  temp¬ 
eratures,  data  from  the  binary  mixture  flows  composed  of  any  two  gases 
from  the  gas  species  He,  Nj  and  CG'2  were  not  included  in  the  statistical 
analysis.  However,  these  data  were  checked  for  fitness  with  the  cor¬ 
relation.  The  r.m.s.  deviation  of  Num  from  the  correlation  line  was  found 
to  be  .0924  when  all  data  below  -  40  were  considered  together. 

5.3  Due  to  the  narrow  range  of  Rem  investigated  for  a  particular 
species  orsped.es  mixture  and,  also  due  to  the  data  scatter,  the  expected 
slight  influence  of  the  small  variation  of  Prandtl  number  (because  of  the 
consideration  of  different  species  and  speciesmixtures  and  also  because  of 
variation  of  temperature)  could  not  be  discei.  ed. 

5.4  Replacement  of  the  usual  temperature  ratio  in  the  temperature 
loading  factor  Ly  a  viscosity  ratio  enabled  a  unique  correlation  to  be 
derived  for  flows  composed  of  gas  species  whose  transport  property  value 
variations  with  temperature  are  different. 

5.5  The  direction  of  the  temperature  loading  effect,  i.e.  in 
the  present  investigation  for  heating  of  cylinders  a  decrease  in  hd/k^ 
with  an  increase  in  temperature  loading,  is  consistent  with  previously 
observed  increase  in  hd/1%  with  an  increase  in  temperature  loading  for 
the  case  of  cooling  of  cylinders. 

5.6  The  Rera  dependency  of  Nu  was  found  to  be  different  from 
that  according  to  King’s  law,  as  also  observed  in  previous  investigations 
with  hot  wires.  However,  due  to  absence  of  previous  data  for  the  case 
of  heating  of  cylinders,  no  direct  comparison  could  be  made. 
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9.?  A  discontinuity  in  the  heat  transfer  curve  was  noted  between 
Re^.  *  hO  and  Re^  *  55.  Unfortunately,  no  data  were  collected  between  these 
values  of  Re^,  and  since  no  specific  aerodynamic  investigation  was  carried 
out  to  determine  a  change  in  the  flow  features  in  this  range,  no  firm 
conclusion  as  to  the  exact  Re^  and  nature  of  the  discontinuity  can  at  present 
be  made. 


5.8  Finally,  a  comparison  of  the  Nujn  values  calculated  for  the 
case  of  heating  of  cylinders  from  correlations  obtained  specifically  for 
cooling  of  cylinders  with  obtained  in  the  present  experiment  (for 
heating  of  cylinders)  clearly  shows  the  need  for  two  separate  correlations 
for  the  two  directions  of  heat  transfer.  The  error  Involved  in  predicting 
heat  transfer  in  one  direction  from  a  correlation  obtained  for  heat  transfer 
in  the  opposite  direction  may  not  be  negligible,  especially  when  significant 
temperature  loading  is  involved. 
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APPENDIX  A 

Expressions  for  the  Calculation  of  Transport  Properties 
for  Pure  Gases  and  Gas  Mixtures 


NOMENCLATURE 


x-1 

x" 


k 


a 

T 

X 

6 

}L 

6 

£l 


SUBSCRIPTS 


nix 


n 

1.2 


Dimensionless  heat  capacity 
Thermal  conductivity 
Monatomic  thermal  conductivity 
Internal  thermal  conductivity 
Boltzman  constant 
Molecular  weight 
Universal  gas  constant 
Absolute  temperature 
Mole  fraction 

Maximum  energy  of  attraction 
Absolute  viscosity 
Zero  energy  collision  diameter 
Reduced  collision  integral 


UNIT 


cal 

cm.  Sec^X 

R 

H 

erg/°X 
g/g  -  mole 
cal/ mole  °X 
°X 


erg. 

g/ca.  sec. 
A 


Components  i  and  j  of  a  mixture 
Entire  mixture 

Number  of  components  in  the  mixture 
Components!  and  2  of  a  mixture 


The  following  expressions  were  used  to  calculate  transport 
property  values  for  non-polar  pure  gases  and  gas  mixtures. 

i)  Viscosity  (Refs,  1.2,3) 
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APPENDD:  A  (Contd) 


or  =  2fc.693  x 
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ii)  Thermal  conductivity  (Refs  1,3) 
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X  '  4-  x" 


K*  =  198.91  x  lu-6 
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X 1  x  0 
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iii)  Viscosity  of  gas  mixtures  (Refs  3,4) 
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V _ /** 
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iv)  Thermal  conductivity  of  gas  mixtures  (Refs  3,4,5) 
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APPENDIX  A  (Con<d) 
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Ths  iim.ensioniess  heat  capacity  i.e.  V'p/R  values  as  functions  of 


temperature  were  cocaine 
(for  viscosity  3rd  therm, 
temperate? 


.V.**  * 


Ref-  6. 


The  reduced  ccllison  integral 


cr.ductivity )  vai  res  at  functions  of  reduced 
d  were  obtained  from  a  table  of  Ref.  ?.  Values  of  the 

r  nr  p* 


rest  of  the  constants  were  taken  from  Table  I  of  Ref.  3 
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FIGURE  2  -  Cooled  Fils  Senaor 
(Photograph  -  Magnification  x  5,  Sketches  -  Not  to  Seale) 
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FIGURE  4  -  Diagram  of  Probe  with  Sensor  Holder  Inserted 
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SLEEVE 
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PROBE  COOLING 

TUBE  BUNDLES 


FIGURE  5  -  Probe  Tip  With  Sensor  in  Operating  Position 
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FIGURE  6  -  Functional  Schematic  of  Circuit 


FIGURE  7  -  Cooled-Film  Circuit  and  Sensor 
Water  Supply  Control  Panels 


(Ts-Tc)°K 


FIG-9  SENSOR  RESISTANCE  Vs  TEMPERATURE  CALIBRATION 
CURVE  OBTAINED  FOR  3  SENSORS. 


FIGURE  10  -  Schematic  of  Experimental  Set-Up 


PIGURE  11  -  Exp«rlMnta\  Arr*ng«Mnt 


FIGURE  12  -  Cloaa-Up  View  of  Plasma  Torch  and  Calibration  Instruments 


FIGURE  13  -  Recording  Station  Arrangement 
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FIG.  14  TYPICAL  INSTRUMENT  SIGNAL  RECORDINGS  FROM  A  SINGLE  RUN 


riGURK  15  ~  Microscope  in  Place  of  Plasma  Torch  for  Instrument  Alignment 
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FLOW  -  N2 
3  T*  -  VARIABLE 
T-  -  CONST 


FLOW  -  He 


T  -  VARIABLE 
T.  -  CONST 


FIG.  16  TABLE  SHOWING  VALUES  OF  CONSTANTS  EVALUATED  BY 
A  LEAST  SQUARE  FIT  OF  THE  DATA  FOR  DIFFERENT 
CONDITIONS  WHEN  A  TEMPERATURE  RATIO  IS  CONSIDERED 
IN  THE  CORRELATION. 


FIG.  17  TABLE  SH0WSN6  VSUje  Or  CONSTANTS  SVAUiVTED  6Y  A 

LEAST  SOUaSE  FIT  OF  THE  DATA  FOB  DFTESEJ.T  CC«€*TCKS 
Wt&i  A  KiNEWUTJC  ViSCCSTTY  RATIO  IS  COJiSQERED  W  THE 
CORRELATION. 
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FIG.I8A 

NUm  V?  R*m  038  A  CIRCULAR  CYLINDER  IN  NITROGEN  CROSS -FLOW  FOR  VARIOUS 
TEMPERATURE  LOftOWGS.  TEMPERATURE  LOADING  EFFECT  OBTAINED  BY  VARYING 
THE  FLO#  TEMPERATURE  AND  KEEPING  THE  CYLINDER  TEMPERATURE  CONSTANT. 
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F1G.I9A 

Num  Vs  Re^  FOR  A  CIRCULAR  CYLINDER  IN  HELIUM  CROSS  -  FLOW  FOR  VARIOUS 
TEMPERATURE  LOADINGS.  TEMPERATURE  LOADING  EFFECT  C3TAINED  SY  VARYING 
T«  AND  KEEPING  Ts  CONSTANT. 
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FIG.  19  B 

Num  (IWiAn)15  Vt  Pem  FOR  HEUUM  CROSS -FLEW  SHOWING  ELIMINATION  OF 
RESIDUAL  TEMPERATURE  LOADING  EFFECT. 
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Nu  Vs  Re  FOR  A  CIRCULAR  CYLINDER  IN  NITROGEN 

m  m  _ 

CROSS -FLOW  FOR  VARIOUS  TEMPERATURE  LOADINGS. 
l*EMPERATURE  LOADING  EFFECT  OBTAINED  BY  VARYING 
THE  CYLINDER  TEMPERATURE  WHILE  KEEPING  THE 


FLOW  TEMPERATURE 

CONSTANT. 
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Num  {I/<D/I/rn)J5  Vs  Rem  rOR  NITROGEN  CROSS -FLOW  SHOWING 
ELIMINATION  OF  RESIDUAL  TEMPERATURE  LOADING  EFFECT. 
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FIG.  2!  a 

Nu  V*  Re  FOR  A  CIRCULAR  CYLINDER  IN  HELIUM  CROSS -FLOW 

m  rfil  t  - - 

FOR  VARIOUS  TEMPERATURE  LOADINGS.  TEMPERATURE  LOADING 
EFFECT  OBTAINED  BY  VARYING  Ts  WHILE  KEEPING  T®  CONSTANT. 


CYLINDER 

DiA  .0152 

cm. 

FLOW 

-  HELIUM 

Ts*K 

too 

Tm 

Ifa 

Vm 

SYMBOL 

387.0 

1556- 

1557 

2044' 

2.C45 

X 

460.0 

1493- 

1.494 

L9I0- 

1920 

o 

532.0 

L436- 

1.437 

1792- 

1802 

a 

FIG.  21 B 


Num '(K./lW5  vs  Rem  FOR  HELIUM  CROSS -FLOW  SHOWING 
ELIMINATION  OF  RESIDUAL  TEMPERATURE  LOADING  EFFECT. 
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!tutnV,  R*m  ***  A  CIRCUUVR  CYUJOER  IN  HELIUM -NITROGEN  MIXTURE 
CROSS-FLOW  FOR  VARIOUS  TEMPERATURE  LOADINGS.  TEMPERATURE  LOACHNG 
EFFECT  OBTAWED  BY  VARYING  T„  AND  KEEPING  T»  CONSTANT 
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FIG.24A 

NumV*  R*m  F0R  A  CIRCULAR  CYLINDER  IN  HEUUM -CARBON  DIOXIDE  MIXTURE 
CROSS -FLOW  FOR  VARIOUS  TEMPERATURE  LOADINGS.  TEMPERATURE  LOADING 
EFFECT  OBTAINED  BY  VARYING  Tm  AND  KEEPING  T»  CONSTANT 
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Num0'./Mt,)-1-  V«  FDR  A  CIRCULAR  CYLINDER  IN  CROSS  -FLOW  OF  VARIOUS  SPECIES  ANO  SPECIE  MIXTURE. 


